Extreme events in two dimensional disordered nonlinear lattices 
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Spatiotemporal complexity is induced in a two dimensional nonlinear disordered lattice through 
the modulational instability of an initially weakly perturbed excitation. In the course of evolution 
we observe the formation of transient as well as persistent localized structures, some of which have 
extreme magnitude. We analyze the statistics of occurrence of these extreme collective events and 
find that the appearance of transient extreme events is more likely in the weakly nonlinear regime. 
We observe a transition in the extreme events recurrence time probability from exponential, in the 
nonlinearity dominated regime, to power law for the disordered one. 
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Nonlinear lattices form a unique laboratory where both 
coherent and incoherent processes appear simultaneously 
and affect their dynamical properties [H- They consti- 
tute prototypical spatiotemporally complex systems that 
can be studied both theoretically and experimentally giv- 
ing a wealth of information on the physics of spatially 
extended complex systems. The additional presence of 
quenched disorder introduces a mechanism for local sym- 
metry breaking that has consequences in the long-time 
dynamics 0] ■ One feature of selforganization of the lat- 
tices relates to the possibility for occurrence of extreme 
events (EE), i.e. the generation of transient or persistent 
structures with large amplitudes that are statistically not 
very significant. In the case of ocean waves, events of this 
type are rogue or freak waves, i.e. waves of considerable 
hight that " appear from nowhere and disappear without 
a trace" Q , usually in relatively calm seas [3, [H ■ Ocean 
rogue waves are attributed to self-focusing effects follow- 
ing the development of Benjamin- Feir instability^. Ex- 
treme waves of this type have also been reported in non- 
linear optical fibers 0, Q , wtical cavities Q , superfiuids 



semiconductor lasers [1^ , 
], soft glass photonic crys- 



[lot , microwave scattering II II 
femtosecond filamentation [13 
tal fiber [ill, ^tc. 

In the present work we focus on two dimensional dis- 
crete nonlinear systems with disorder that have been re- 
alized experimentally [l^ and address the question of 



the conditions for generation of very large, i.e. extreme, 
discrete wave events in the process of spatiotemporal evo- 
lution. While the production of EEs is mediated by mod- 
ulational instability (MI) the subsequent evolution shows 
a complex dynamics that leads to generation as well as 
destruction of waves of large magnitude. This dynamics 
is directly induced by the interplay of nonlinearity and 
disorder; we are seeking the regimes where EE appear- 
ance is the most frequent as well the statistics of their 
appearance. In recent work we found that integrability 
properties of the lattice do play a role in the generation of 
EEs [l3|- The present work, focusses on the issue of dis- 
crete lattice EEs from a broader perspective and through 
a physically realizable model, viz. that of the Ander- 
son disordered Discrete Nonlinear Schrodingcr (DNLS) 
equation [l^ [l^ . We address two experimentally rele- 
vant questions, viz. what is the optimal regime for EE 
generation and what is the recurrence time statistics of 
these events once generated. From this analysis a clearer 
picture of the spatiotemporal complexity of the system 
emerges. 

We consider the dynamics of a strongly coupled excita- 
tion in a two dimensional tetragonal lattice with diagonal 
disorder, or, equivalently, wave propagation in a two di- 
mensional array of evanescently coupled optical nonlin- 
ear fibers with random index variation, both described 
through the disordered DNLS equation, viz. 



dt 



+ J{'4^n+l,m + V'n — l,m H" V'n.rn+l ~l~ Wn.m—l ) + 7l'0 
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where n,m = 1,...,N, '>pn,m is a probability (or wave) 
amplitude at site (n, m) , J > is the inter-site coupling 
constant accounting for tunnelling between adjacent sites 
of the lattice (corr. evanescent coupling), 7 > is the 



nonlinearity parameter that stems from strong electron- 
phonon coupling (corr. Kerr nonlinearity), while en,m, 
is the local site energy (related to the fiber refractive 
index), chosen randomly from a uniform, zero-mean dis- 
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tribution in the interval [-W/2,+W/2]. Equation ^ 
serves as a paradigmatic model for a wide class of phys- 
ical problems where both disorder and nonlinearity are 
present. For 7 0, Eq.dT]) reduces to the 2D Ander- 
son model while in the absence of disorder (e„^„i — 0), 
it reduces to the DNLS equation in 2D that is gener- 
ally non-integrable. Equation ([!]) conserves the norm 

Pn = J2n=i Sm=i l^n.mP' ^nd the Hamiltoniau H, cor- 
responding to total probability (corr. input power) and 
the energy of the system, respectively. Wc will con- 
sider the evolution of an initially uniform state that is 
slightly modulated by periodic perturbations in order to 
facilitate the development of MI [13], [l^. Specifically, 
we use the initial state 4'n,m = ("00 + ^'0n,m) exp(— z/it), 
where ipo is a constant amplitude of the continuous wave 
(CW) (averaged over disorder, with (e„,„i) = 0), and 

n^m = OLn,m + ibn,m En'e Small complcx pcrturbatious 
that modulate the constant CW solution. The MI in- 
duces nonlinearly localized modes that are, however, 
modified by the presence of the quenched disorder. In 
the time scale of the numerical study, we note that the 
presence of disorder induces some form of additional en- 
ergy redistribution among the sites of nonlinear lattice. 

We adopt the criterion employed in oceanography and 
define as an extreme wave one that has height Hext (here 
the positive difference of minimum-to-maximum ampli- 
tude iV'n.ml in different time steps) with Hext ^ 2.2_ffs, 
where iJ^, is the significant wave height, i.e. average 
height of the one-third highest waves in a given sample 
0. We integrate Eqs. ([!]) with periodic boundary con- 
ditions using a 6— th order Runge-Kutta [l3| solver for 
several different values of W and 7 (J = 1) using typ- 
ically a lattice with = 41 (we have checked results 
for = 81 as well) and periodic boundary conditions. 
During relatively long time (typically we used maximum 
time of 1 X 10'^ units or equivalently approximately 500 
coupling lengths) the system selforganizes and wc ob- 
serve localized structures appearing on different sites sur- 
rounded by irregular, low-amplitude background. Some 
of these structures are in the form of breathers, either 
pinned or mobile, while some others are transient. The 
complete amplitude statistics for the observed time in- 
terval is shown in Fig. ([1^) for different disorder and 
nonlinearity parameters; we observe a distinct, Ray leigh- 
like distribution that vary depending on the parameters. 
Any of the states, longer lived or transient, that appear in 
tail of the distribution are candidates for extreme events. 

In order to quantify the onset of extreme events in the 
lattice we use several measures, viz. the inverse partici- 
pation ratio, the probability for extreme event occurrence 
as well as first appearance and recurrence EE times. The 
inverse participation ratio 
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FIG. 1: (color online) Probability amplitude distribution 
Ps as a function of amplitude \ip\ for different nonlinearity 
strengths and W = 2 (a). The curves have a Rayleigh-like 
behavior, (b) Effective average width of excitations formed 
in the lattice as a function of time for different disorder levels 
for 7 = 0, 2. We note that for fixed nonlinearity, increase in 
the amount of disorder induces delocalization. 



determines the effective confinement of initially extended 
excitations. We define i^eff = P^^^; this effective width 
gives the average spatial extent of the structures gener- 
ated after the MI (Fig. (Pd)). 

Both disorder and nonlinearity, when each acts alone, 
favor wave localization in the lattice. When they are 
present simultaneously, quenched disorder dominates the 
early dynamics since MI develops more slowly, at least for 
relatively small 7's. In this regime, where both nonlinear- 
ity and disorder are present, the Anderson-like localized 
states decay spatially in a fashion that still permits local 
energy redistribution, while at relatively longer times a 
more stable state is reached. We note from Fig. ([Ud) 
that ujef f saturates to lower values for non-zero nonlin- 
earity compared to the value in the corresponding lin- 
ear disordered system. This tendency is compatible with 
the findings of Schwartz et al. who observe experimen- 
tally that increased self-focusing enhances localization 
(3]. On the other hand, we also see in the same fig- 
ure that increasing disorder for fixed nonlinearity results 
to larger We//, leading thus to reduction of localization. 
This opposite tendency of delocalization induced by dis- 
order in the presence of nonlinearity can be attributed 
to the partial destruction of the pinned, highly localized 
nonlinear states by the presence of disorder and the as- 
sociated weak redistribution of energy. 

In order to find the regimes that favor EE production 
we calculate numerically Pee i the probability for EE gen- 
eration with h > Hext = 2.2Hs; it is shown in Fig. [5] 
as a function of 7 for three different levels of disorder. 
From Fig. [2] we observe that increasing the strength of 
nonlinearity results in an increase in the EE probability 
Pee, reaching generally a maximum, while at high non- 
linearity parameter values ( 7 > 2 ), the EE probability 
reduces for all observed disorder levels. We note that 
for small disorder the curve has additionally a secondary 
maximum while the occurrence of the primary maximum 
value shifts to large 7-values for larger amounts of disor- 
der W. Further increase of nonlinearity leads to the de- 
cay of the extreme event generation probability, however, 
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with a tail that is substantially higher for the larger disor- 
der case. For vanishingly small nonlinearity, Pee is small 
but with appreciable values; we find Pge = 0.47, 0.44 for 
disorder levels with W = 2, 4, respectively. These re- 
sults indicate that the most favorable conditions for the 
EE generation during propagation through in the 2D lat- 
tice systems are relatively high level of disorder compared 
to the linearized bandwidth and presence of nonlinearity. 
The interplay of disorder and nonlinearity results in one 
or more maxima in the probability of EE generation for 
certain values of the nonlinearity strength. The first max- 
imum which is the most pronounced one for relatively 
small disorder and weak nonlinearity can be correlated 
with high EEs probability in the near integrable system 



pointed out in the ID context [la]. On the other hand, 
for arbitrary disorder the localizing effect of nonlinearity 
seems to be responsible for the appearance of the second 
maximum at 7 « 2 with a value Pee ~ 0.7%. 




FIG. 2: (color online) Extreme event height probability Pee 
as a function of the nonlinearity strength 7 for several lev- 
els of disorder. Arrows show positions of local probability 
maximums. 

In order to probe deeper into the genesis and appear- 
ance statistics of EEs in discrete, disordered nonlinear 
lattices we focus on the probability distribution P^ of 
the recurrence time r of the extreme events [20|, 211. We 



study two quantities, the mean recurrence time R of an 
EE as well as the mean time of first appearance (TFA). 
As these times depend on the specific value of the pa- 
rameter q, the wave height threshold, we first investigate 
their dependence on q. Numerics show that the slope of 
the curves of the recurrence time P as a function of the 
threshold q is smaller when only disorder is present com- 
pared to the case of simultaneous presence of disorder 
and nonlinearity for all levels of disorder (Fig. [3^). The 
mean recurrence time R increases with q for all values of 
W and 7 (Fig. [3]). The increase is faster for lower disor- 
der level (Fig. ^jji.) and stronger nonlinearity (Fig. (SJd). 



This dependence in the strong nonlinearity-weak disor- 
der regime is related to the favorable conditions for the 
creation of highly pinned, immobile localized structures, 
a fact that increases dramatically the mean recurrence 
time. Qualitatively similar behavior we found also for the 
TFA-q dependence. These studies show the complexity 
of the recurrence as well as first appearance phenomena 
and their direct manifestation in time average quantities. 
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FIG. 3: (color online) The average return time R of extreme 
events as a function of the threshold amplitude q for (a) var- 
ious disorder levels with 7 = (circles) and 7 = 2, (b) fixed 
level of disorder {W = 2) and several nonlinearity strengths 
(including zero nonlinearity). 

We calculate P^ , the return time probability for differ- 
ent disorder levels and nonlinearity values in the following 
way. For a given threshold value q {q > Hext) wc scan 
the lattice and find an event at a given location with am- 
plitude larger than q; we register as recurrence time the 
time between this event and a subsequent above thresh- 
old one that appears in the same lattice location . We 
follow this procedure repeatedly up to maximum time 
and construct histograms for different parameter values, 
all scaled by the average return time R in each parame- 
ter regime; the outcomes for the various non normalized 
return probabilities are shown in Fig. |4]for several values 
of disorder level and nonlinearity for a given g- value. In 
the purely linear but disordered regime we find that P^ 
is a power-law function, viz. Pr = [a + b{r/ R)]~^ , with 
/? w 1.34 q = 0.66; while the parameters depend on the 
threshold value q, the functional form is stable. The addi- 
tion of nonlinearity to a disordered lattice reduces some- 
how the possibility of occurrence of EEs. Remarkably, 
the Pr vs. r/R dependence of the curves changes gradu- 
ally from power-law, in the linear, disordered regime, to a 
double exponential at intermediate nonlinearity values to 
single exponential for a pure nonlinear lattice with no or 
very little disorder. Furthermore, the decay of the curves 
is more rapid for higher nonlinearity strength. The ob- 
served transition from a power law, in the disordered 
dominated regime, to exponential, in the nonlinearity 
dominated one, is linked to the behavior of the tail of 
the amplitude probability distributions. 

We relate our relatively short-time results to expected 
long-time wave-packet spreading regimes summarized in 
Ref. [2^ and [2^ (one-dimensional DNLS with disorder 
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FIG. 4: (color online) Return time probability Pr (not nor- 
malized) as a function of r/R, the return time scaled by the 
corresponding average return time, for the purely disordered 
case (power law fit - black curve) , purely nonlinear case (expo- 
nential fit - red line) and intermediate cases (double exponen- 
tial fits). We use 5 = 0.66. In box the lines which designate 
the first two cases are shown separately. 
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FIG. 5: (color online) Different regimes of the spreading of 
wave-packets in the eflfective parameter space {W, S) [2^. [23|. 
Lines represent regime boundaries 5 k, d and 5 = 2. Empty 
squares denote parameters for which the local maximums of 
Pee are found 



and its generalization to the higher dimensional prob- 
lem, respectively). Different regimes are obtained for 
(i) 5 > 2, onset of selftrapping, (ii) d < 5 < 2, strong 
chaos and (iii) 5 < d for the effective quantities d and 
5. The latter are the average frequency spacing of the 
nonlinear modes within a localization volume given by 
d K, A/F = (8 J -I- W)lijJeff while the nonlinear fre- 
quency shift 5 is taken to be (5 ~ 7. The selection of 
regimes in p^[23j was done taking into account the inten- 
sity of interaction among the nonlinear modes; the latter 



increases with the nonlinearity strength up to the high 
nonlinearity (here 5 « 2) when the strong self-trapping 
results in a creation of isolated strongly pinned high am- 
plitude breathers. The comparison is approximate but 
leads to interesting observations (see Fig. [5]). The first 
local maximum in the Pee found for small nonlinearity 
strength is located in the weak chaos regime. Its max- 
imal value Pee ~ 0.8% is observed for small W . The 
second, broader maximum of Pee for all disorder cases 
is located in the strong chaos regime relatively close to 
the border lines with neighboring regimes (Fig. On 
the other hand, we may associate the power-law decay 
of the Pj. to the weak chaos regime, and the exponen- 
tially like one with the self-trapping regime. Therefore, 
transient EEs are more probable in the regime of weak 
chaos, while the long-lived EEs (high amplitude strongly 
pinned breather structures) are dominant structures in 
the regimes of strong chaos and self-trapping. This as- 
sociation enables us to relate the first local maximum 
in Pee to the weak interaction of transient EEs induced 
by disorder while the second, broader maximum, to the 
appearance of longer-lived breathers resulting from the 
energy redistribution through the strong interaction be- 
tween nonlinear modes. 

We have calculated several statistical measures for 
wave propagation in disordered and nonlinear lattices de- 
scribed by the DNLS in two-dimensions with respect to 
EEs generation. The interplay of two distinct mecha- 
nisms of energy localization leads to several unexpected 
results. The calculations indicate that EEs may be gen- 
erated with appreciable probability also in linear disor- 
dered lattices (but not in linear non-disordered lattices) 
for the chosen type of initial excitations. 

When the relative strengths of nonlinearity and dis- 
order vary, the average return probability Pr of EEs 
as a function of r/R changes between the two limiting 
regimes, i.e., disorder without nonlinearity and nonlin- 
earity without disorder corresponding respectively to a 
transition from power-law to exponential . These findings 
are illustrated in Fig. |3]for certain parameter values. The 
change occurs gradually between the two regimes passing 
through a sequence of states that have more complex re- 
turn time behavior. Moreover, the calculated probability 
for EE generation indicates that there are two parameter 
regimes of the disorder level and nonlinearity strength 
which favor EE generation, as shown in Fig. ([2]). We 
found that a relatively low level of disorder (compared 
to the bandwidth) leads to maxima in the probability for 
EE production in the relatively weak nonlinearity regime. 
The second, broader, maximum found for higher non- 
linearity does not depend significantly on the disorder 
level. For very high 7's the probability of EE generation, 
although small, depends strongly on the level of the dis- 
order. The fact that EE generation is favorable in the 
weakly nonlinear regime appears to be compatible with 
previous findings in discrete lattices, where the highest 
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probability for E E g eneration is observed in the nearly 
integrable limit [l6|. On the other hand, the appear- 
ance of the second parameter regime which favors the 
EE generation is correlated with the localization mecha- 
nism governed by nonlinearity |2J| . We also note that the 
mechanism responsible for the probability for EE gener- 
ation depends additionally on the specific initial condi- 
tions used, as seen also in other approaches [2^. 
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